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It is known that in a hydrodynamic Taylor-Couette system uniform rotation or a rota¬ 
tion law with positive shear (‘super-rotation’) are linearly stable. It is also known that 
a conducting fluid under the presence of a sufficiently strong axial electric-current be¬ 
comes unstable against nonaxisymmetric disturbances. It is thus suggestive that a cylin- 
dric pinch formed by a homogeneous axial electric-current is stabilized by rotation laws 
with d D. j d R > 0. However, for magnetic Prandtl number Pm f 1 and for slow rota¬ 
tion also rigid rotation and super-rotation support the instability by lowering their critical 
Hartmann numbers. For super-rotation in narrow gaps and for modest rotation rates this 
double-diffusive instability even exists for toroidal magnetic fields with rather arbitrary 
radial profiles, the current-free profile B n> oc 1 /R included. - For rigid rotation and for 
super-rotation the sign of the azimuthal drift of the nonaxisymmetric hydromagnetic insta¬ 
bility pattern strongly depends on the magnetic Prandtl number. The pattern counterrotates 
with the flow for Pm <C 1 and it corotates for Pm 1 while for rotation laws with 
negative shear the instability pattern migrates in the direction of the basic rotation for all 
Pm. 

An axial electric-current of minimal 3.6 kAmp flowing inside or outside the inner cylin¬ 
der suffices to realize the double-diffusive instability for super-rotation in experiments us¬ 
ing liquid sodium as the conducting fluid between the rotating cylinders. The limit is 11 
kAmp if a gallium alloy is used. 
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I. INTRODUCTION 


A well-known instability of toroidal fields is the magneto hydrodynamical pinch-type current- 
driven instability which is basically nonaxisymmetricA. The toroidal field becomes unstable if a 
certain magnetic field amplitude is exceeded depending on the radial profile of the field which 
forms the electric-current pattern. It is also known that for unity magnetic Prandtl number a global 
rotation of the system increases the critical field amplitude. The latter is strongly reduced, however, 
if the rotation decreases outwards (i.e. dJ?/d R < 0, ‘sub-rotation’). The formal reason is that sub¬ 
rotation becomes (Rayleigh-) unstable even in the hydrodynamic regime if it is steep enough. More 
important is the existence of a nonaxisymmetric instability for such rotation laws even for current- 
free toroidal fields oc 1 /R) which we have called azimuthal magnetorotational instability 
(AMRI 2 A), it appears for all values of the magnetic Prandtl number for rather low Hartmann 
numbers but for large magnetic Reynolds numbers of the basic rotation. This phenomenon also 
explains the general destabilization of toroidal fields by rotation laws with 1? decreasing outwards. 

The question arises about the role of ‘super-rotation’, i.e. rotation laws with di?/d R > 0, 
which are linearly stable in the hydrodynamic regime^. The nonlinear behavior is less clear as 
some Taylor-Couette experiments have shown instability in this regime-^. Superrotation cannot 
be destabilized by the standard magnetorotational instability with axial external fields. Inspired by 
the discovery of the axisymmetric helical MRI a WKB method for inviscid fluids in current-free 
helical background fields has been applied providing two limits of instability in terms of the shear 
in the rotation law- An upper threshold suggests a magnetic destabilization of super-rotating flows 
for very strong positive shear. A similar phenomenon has been reported by Bonanno & Urpin 
(2008) resulting from a local analysis for a helical field under the influence of super-rotation^. 
Later it has been shown with a dispersion relation for inductionless fluids (see the Appendix) that 
the stability curve does not cross the line representing the differentially rotating pinch formed by 
uniform electric-current suggesting instability for both signs of shearAA. 

By means of a corresponding approximation Acheson (1978) showed that for fast-rotation the 
current-driven instability of toroidal fields may be stabilized by positive shear 12 -. If this is true 
we expect in the solar low latitudes where in the bulk of the convection zone the equatorial i? 
increases outwards that the toroidal field is stabilized and can be amplified to much higher values 
than it would be true for the opposite rotation law. Contrary to that a rotation law with negative 
shear - as it exists in higher solar latitudes - strongly destabilizes the fields so that they cannot 
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reach high amplitudes. It is shown here by use of a simplifying cylinder geometry that indeed for 
not too small magnetic Prandtl numbers super-rotation stabilizes toroidal magnetic fields while 
sub-rotation strongly destabilizes toroidal magnetic fields. On the other hand, small magnetic 
Prandtl number and slow rotation of any rotation law - including rigid rotation - lead to lower 
critical magnetic field strengths than needed for destabilization at i? = 0. We shall show in the 
present paper that for di?/d R > 0 and for slow rotation the relaxation of the excitation conditions 
compared with the resting container belongs to the double-diffusive phenomena which disappear 
if the the molecular viscosity equals the molecular resistivity. 

A Taylor-Couette container is considered which confines a toroidal magnetic field with ampli¬ 
tudes fixed at the cylinders which may rotate with different rotation rates. The gap between the 
cylinders is considered as variable. Normalized with the outer radius R oat the inner radius R m is 
> 0.5. The cylinders are unbounded in axial direction. 

The fluid between the cylinders is assumed to be incompressible and dissipative with the 
kinematic viscosity v and the magnetic diffusivity //. Derived from the conservation of angular 
momentum the rotation law Q(R) in the fluid is 

f?(f?) = a + 4 (!) 


with 

a = f?in, b = J—lLRln^, (2) 
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fi m and i? out are the imposed rotation rates of the inner and outer cylinders. After the Rayleigh 
stability criterion the flow is hydrodynamically stable for // > rf n . We are only interested in 
hydrodynamically stable regimes so that // > rf n should always be fulfilled. Rotation laws with 
di?/d R > 0 are described by fi > 1 while rotation laws with di?/d R < 0 are described by 
/i < 1. Rigid rotation means n — 1. Hydrodynamical flows with rigid rotation or super-rotation 
are always linearly stable. 

Also the possible magnetic profiles are restricted. The solution of the stationary induction 
equation without flows reads 


— AR 


( 4 ) 
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(in cylinder geometry, see Roberts 1956^, Tayler 1957—) where A corresponds to a uniform axial 
current everywhere within R < R 0 ut . The quantity B out /B m measures the variation of B 0 across 
the gap. For fields after © it is simply B out = P in /r in . If the axial electric-current only exists 
inside the inner cylinder then the solution of the induction equation instead of © is 

B<j, oc (5) 

In the present paper the stability characteristics of the MHD system are due to the instability of 
the field © under the influence of super-rotation. In order to compare the results for the standard 
profile © with those for the field which is current-free in the fluid, the profile © has only been 
used below for data given in Fig. |4j We know that the nonaxisymmetric Tayler instability (TI) also 
exists for resting fluids with threshold values which do not depend the magnetic Prandtl number 

Pm — (6) 

V 

the value of which, however, has an essential influence on the excitation of the TI under the in¬ 
fluence of rotation^. For fast rotation, a narrow gap, Pm ^ 1 and strong shear we shall present 
instability maps which hardly differ for various magnetic profiles, the vacuum fields of AMRI 
included. For these solutions, therefore, the importance of the electric-current inside the fluid dis¬ 
appears and the instability gets its entire energy from the differential rotation. This conclusion will 
be supported by the inspection of the associated wave numbers and drift velocities of the various 
nonaxisymmetric instability patterns. 

Both AMRI and (resting) TI have recently been realized in the MHD laboratory using the 
liquid eutectic alloy GalnSn with Pm = 1.4 • 10 -6 as the conducting fluid^^. If the results shall 
be applied to turbulent media like the stellar convection zones then the magnetic Prandtl number 
must be replaced by its turbulence-induced values which are much larger—. In the upper part of 
the solar core the molecular value is about Pm ~ 0.065^. 

II. EQUATIONS 

The dimensionless incompressible MHD equations are 

f c)u \ 

Re f — + (u ■ V)w j = —VP + A u + Ha 2 curl B x B, 

dB 

Rm— = curl(« x B) + A B, (7) 
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with div u = div B = 0 and with the Hartmann number 


Ha = 


B-, n D 


( 8 ) 


yJ\M>pVV) 

D = yj -Ri n (-R ou t — R m ) is used as the unit of length, ft-- 1 as the unit of the time, r)/D as the unit 
of velocity and B m as the unit of magnetic fields. In this notation the angular velocity of the global 
rotation ft at the inner cylinder equals Rm. The Reynolds number Re is defined as 

fLD 2 


Re = 


v 


( 9 ) 


and the magnetic Reynolds number as Rm = Pm Re. It is also useful to work with the mixed 
Reynolds number 

Rm = \/ReRm (10) 

which is symmetric in v and ?/ as it is the Hartmann number. Its ratio to © is called the magnetic 
Mach number Mm which measures the rotation rate in comparison with the Alfven frequency 

ft A = B in /y/n 0 pD 2 , 

ft\ i 


Mm = 


Rm 


(ID 


J?a Ha 

We always use no-slip boundary conditions for the velocity u R = = u z = 0. The material 

of the cylinders is assumed as made from perfect conductors, or in some other cases made from 
perfect insulators. For the conducting walls the fluctuations b have thus to fulfill the conditions 
d&^/df? + b^/R = b R = 0 at both R in and R out . Mathematical details about the much more 
complicated vacuum boundary conditions and the used numerical codes can be found in previous 
publications. The time-tested code for the linearized equations^ solves the eigenvalue problem for 
the Fourier modes exp(i(cct + kz + mo )) where k is the axial wave number and m the azimuthal 
mode number. 

The nonlinear simulations have been done with our reliable time-stepping code. It works with 
an expansion of the solution in Fourier modes in the azimuthal direction generating a sample of 
meridional problems each of which is solved using a Legendre spectral element method 21,22 . 


III. A DOUBLE-DIFFUSIVE INSTABILITY 


We start with solid-body rotation of a container. It is known that this flow with © belongs to the 
class of magnetized flows where the radial profiles of the global velocity and the global magnetic 
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field are identical. Chandrasekhar (1956) has shown that all flows of this class are stable™ One 
can also show for the nonideal MHD flows that all lines of marginal instability for very small Pm 
are identical in the Ha-Re planed. They are given for various magnetic Prandtl numbers in Fig. 
Q] (left panel). For small Pm the lines for Pm < 10 cannot be separated optically. It is also 
demonstrated that the critical Hartmann number Ha Tay for Re = 0 does not depend on Pm. The 
curves, however, for Re > 0 behave different. For Pm < 1 they turn to the left while for Pm > 1 
they are turning to the right. In the former case the instability is supported by the rotation and in 
the latter case it is (strongly) suppressed. Note that the lowering of the Hartmann number by global 
rotation only exists for Pm 1 and for slow rotation with Mm <C 1. For those parameters “the 
stabilizing effect of global rotation is greatly reduced”, as it already has been formulated for this 
kind of double-diffusive problems™ A very similar behavior also appeared for the excitation of 
axisymmetric modes for helical background fields ( B z B (t) ^ 0) under the influence of differential 
rotation which are also stable for u — q — 0 and their critical eigenvalues are lowered for Pm ^ 1 
compared with those for equal diffusivities™ 

For faster rotation all curves are turning to the right so that the statement finally becomes correct 
that generally the rotation suppresses the Tayler instability. 
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FIG. 1. Instability of the m = ±1 mode for rigid rotation. The numbers represent the magnetic Prandtl 
number Pm. Left panel: The lowering of the critical Hartmann numbers by the global rotation for Pm / 1 
Note the maximal relaxation of the Ha existing for Pm = 0( 1CP 2 ) but not for Pm —>• 0. Right panel: The 
drift rates oj dr /Q of the modes for marginal instability for small Pm and large Pm. The axial wave numbers 
k for all points at the lines do hardly vary. r- m = 0.5, perfect-conducting boundaries. 


Interesting is also the behavior of the drift rate cu dT as the real part of the frequency oo of the 
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Fourier mode of the instability pattern, normalized with the rotation rate of the cylinders. It can be 
positive or negative. Because of the definition 


CUdrf^out 


m 


( 12 ) 


the azimuthal migration of the instability pattern has the opposite sign of cudr- The solutions for 
the modes with m — 1 and m — — 1 have the same eigenvalues Re and Ha but the opposite 
signs of cudr- After ( ITU) they have thus the same frequency of migration in ^-direction. Figure Q] 
(right panel) shows that the drift rate hardly depends on the Hartmann number but it is strongly 
directed by the magnetic Prandtl number Pm. It is positive for small Pm which with (fl2l) leads to 
an azimuthal migration of the instability pattern opposite to the cylinder rotation but for Pm > 1 
it changes the sign so that the instability pattern rotates in the same direction as the cylinders do. 
The relaxation of the critical Hartmann number for Pm ^ 1 is not indicated by the drift rates. 

Figure [2] summarizes the influence of nonuniform rotation on the excitation of the TI against 
nonaxisymmetric modes with m = 1. It is r in = 0.8 and the cylinders are made from insulating or 
perfect-conducting materials. The critical Hartmann number for resting cylinders is Haxay = 250 
for vacuum boundary conditions and HaT ay = 290 for perfect-conductor conditions. Again, these 
values do not depend on the magnetic Prandtl number as here demonstrated here for the two 
examples with Pm = 1 (left panel) and Pm = 10~ 5 (right panel). One finds that the TI can be 
excited more easily for vacuum boundary conditions than for perfect-conducting cylinders. 

The dotted lines in Fig. [2] are the lines of marginal instability for the sub-rotation law with 
/a = 0.5. In the narrow gap and for fast enough rotation such a rotation law is linearly unstable 
without magnetic field. For not too small Pm the magnetic field even destabilizes such a steep 
sub-rotation law so that for finite Hartmann number the critical Reynolds number is always lower 
than 160 which is the critical value for Ha = 0. The influences of the boundary conditions and the 
magnetic Prandtl number on the stability/instability of sub-rotation laws is only small. 

The solid lines are due to rotation laws with positive radial shear. Their behavior strongly 
depends on the value of the magnetic Prandtl number. For Pm = 1 (left panel) super-rotation 
(/i = 4) acts stabilizing (dRe/dHa > 0 everywhere) while for small magnetic Prandtl number 
(Pm = 10 -5 , right panel) it acts destabilizing (dRe/dHa < 0 in the lower part of the diagram). 
Under the presence of rotation with positive shear the electric-current becomes unstable for lower 
Hartmann numbers than for the resting pinch. This ‘subcritical excitation of the TI for super¬ 
rotation is insofar interesting as TC-flows with positive shear are prominent examples of stable 
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FIG. 2. Instability map for the m = ±1 mode for Pm = 1 (left) and Pm = 10~ 5 (right). There are 
examples for sub-rotation (// = 0.5, dotted lines), rigid rotation (q = 1, dashed lines) and super-rotation 
(q = 4, solid lines). The lowering of the critical Hartmann numbers for super-rotation only appears for 
Pm / 1. The boundary conditions are those for vacuum or for perfect conductors, r m = 0.8. 


hydrodynamic flows^. Such a very stable configuration can even be destabilized by a magnetic 
field which is weaker than the critical field for the TI. This phenomenon only exists for slow 
rotation since for fast rotation all nonaxisymmetric magnetic instabilities are suppressed by any 
sort of differential rotation. The magnetic Mach number Mm for the subcritical excitation by 
super-rotation in the right panel of Fig. [2] is (only) of order 10~ 3 . Note that a rotation profile 
with q > 1 needs much higher Hartmann numbers to be destabilized than a rotation law with 
q < 1 and sufficiently high enough Reynolds number. With other words, the toroidal field which 
is induced by a super-rotation can become much stronger than the toroidal field which is induced 
by a sub-rotating J?-profile. This basic finding should have implications for the electrodynamics 
of rotating stars. With other words, the instability requires rotation laws with negative shear in 
order to exist for large magnetic Mach numbers. 

The dashed lines describe the influence of rigid rotation which for fast rotation always acts 
stabilizing, i.e. dRe/dHa > 0. This effect, however, is much weaker for small Pm than for 
Pm = 1. In the latter case the instability is suppressed for all slow rotation rates. The stabilization 
depends on the value of the magnetic Prandtl number, it is strongest for Pm = 1. 

In summary, we have found that for Pm = 1 rigid rotation and super-rotation suppresses the TI 
for all Re. Rigid rotation and super-rotation support the TI for Pm ^ 1 but only for slow rotation 
by lowering the Hartmann number, Ha < IhiTay, below the value which holds for the resting 
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pinch. The lowering only exists if the two molecular diffusivities have different values - no matter 
which is larger or smaller - typical for a double diffusive instability^’ 2 -. 

From Fig. [2] one also finds that the boundary conditions do not play a minor role. For perfect¬ 
conducting cylinders the lowering of the Hartmann numbers is much stronger than than for cylin¬ 
ders of insulating material. In order to study the influence of the geometry on the instability we 
shall consider in more detail the two TC-flows with a rather narrow gap (r in = 0.95) and with a 
wide gap (r in = 0.5) . 

IV. NARROW GAP 

For a narrow gap the influence of the magnetic Prandtl number on the Tayler instability in 
a container with various rotation profiles shall be studied. All the considered rotation laws are 
hydrodynamically stable. For a gap with r in = 0.95 Fig. [3]gives the results for Pm = 0.1, Pm = 1 
and Pm = 10. The critical Hartmann number for TI without rotation is HaT ay = 3060. For 
Pm = 1 rigid-body rotation and super-rotation of any Reynolds number are always stabilizing, i.e. 
Ha > 3060. Only sub-rotation leads to Ha < 3060. The differences for both sub-rotation and for 
super-rotation here only appear for rather low Reynolds numbers. 

It is also worth to mention that the lines of marginal instability for rigid rotation and for super¬ 
rotation always lie below the line Ha = Rm, i.e. even slow rotation stabilizes the TI for Pm = 1. 
For fast rotation the TI only exists under the presence of differential rotation with negative shear. 





FIG. 3. Stabilization and destabilization by super-rotation in a narrow gap for Pm = 10 (left). Pm = 1 
(middle) and Pm = 0.1 (right). The modified Reynolds numbers (ITOl ) are used for an easy definition of 
the unity magnetic Mach number Mm (dashed lines). For comparison one example for sub-rotation is 
given which itself is hydrodynamically stable. The curves arc marked with their value of //. r m = 0.95, 
perfect-conducting boundaries, m = ±1. 
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A. Subcritical excitations 


To discuss the results for Pm ^ 1 it makes sense to use the modified Reynolds number Rm for 
the characterization of the basic rotation. For both Pm > 1 (Fig. [2 left panel) and Pm < 1 (Fig. [3] 
right panel) also the rotation laws with positive shear lead to subcritical excitations if the rotation 
is slow enough. The magnetic Mach number which measures the rotation of the (inner) cylinder 
to the Alfven frequency for the subcritical excitation and for both magnetic Prandtl numbers is 
Mm ~ 0.05. Again the curves for rigid rotation and for super-rotation are located below the line 
Mm = 1. Again, for Mm > 1 the TI needs the action of a sub-rotation law with negative shear. 

For sufficiently fast rotation the super-rotation laws are always stabilizing. The super-rotation 
for small magnetic Prandtl numbers is much more stabilizing than that for high magnetic Prandtl 
numbers. For Pm = 10 the stabilization by super-rotation is even weaker than that of rigid 
rotation. It is often the rule for magnetic instabilities that large Pm destabilize nonuniform rotation 
while small Pm stabilize the flows. The formal reason for this phenomenon can be realized in 
Fig. [2] where the bifurcation curve for super-rotation for small Pm moves to the left of the line for 
rigid rotation rather than to the right as for Pm = 1. 

It is known that rotation laws with negative shear (here // = 0.92) behave strongly destabilizing 
the flow. The domain of stability in Fig. [3] is again larger for small Pm. For sufficiently fast 
rotation also the lines of marginal instability for sub-rotation turn to the right stabilizing the system 
as rotation laws of strong shear of both signs do always erode nonaxisymmetric magnetic patterns. 

The question arises about the possible existence of a minimum Hartmann number for steeper 
and steeper super-rotation laws. The existence of such a limit is suggested by the suppression of 
a nonaxisymmetric magnetic field by differential rotation which should grow with growing values 
of shear. The line of marginal instability can never touch the vertical axis as without magnetic 
field super-rotation always behaves stable. Figure 0] shows very close lines for // = 4, // = <8 
and even fi = 128 so that the minimum Hartmann number HaMm can be estimated with strongest 
super-rotation as smaller by a factor of three compared with HaT ay = 3060. For the very small 
magnetic Prandtl number used for Fig. 0] (left) the numerical value of HaMin/HaTay is astonishing 
small. For large Pm (right panel of Fig. 01 Pm = 10) the subcritical excitation also occurs with 
HaMin/HaTay even smaller. For larger Reynolds number almost all curves (except the curve for 
rigid rotation) are identical, they depend on numerical values of shear and current only in a very 
little manner. Compared with the curves for very small Pm, however, the curves have a different 
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FIG. 4. Instability maps for super-rotation for small Pm (Pm = 10 -5 , left) and large Pm (Pm = 10, 
right), the lines marked with their values of /i. The curves are given for the two extreme radial profiles © 
and © by solid and by dashed lines. Here, the Reynolds numbers are defined with the outer rotation rate. 
n n = 0.95, perfect-conducting cylinders. 

There is another striking feature plotted in Fig. |4l In the domain where the lines are almost 
vertical the dependence of the critical Reynolds number on the critical Hartmann number is ex¬ 
tremely weak. It is shown in this plot that there even the dependence of the curves on the radial 
profile of B = B rh (R) is weak. The dashed lines in Fig. [4]represent the instability for the field 
(0) which is current-free between the cylinders within the fluid. These curves, therefore, can never 
cross the horizontal axis where Re = 0. For fast rotation, however, they almost coincidence with 
the solid lines for the marginal instability of the flow with axial current. Surprisingly, for strong 
shear and fast rotation the presence of the electric-current becomes irrelevant for the occurrence of 
instability. One can show that all possible radial profiles of B a between © and ([5]) provide more 
or less the same instability curves in this domain of the bifurcation map revealing that the differ¬ 
ential rotation for Pm ^ 1 is able to deliver the entire energy for the maintenance of the instability 
patterns and the magnetic field only acts as a catalyst. This phenomenon is already known from 
AMRI for sub-rotation but here, for super-rotation, it only works for v j- tj. 

The close relatedness of both the instabilities for the lowest Hartmann numbers is obvious. It 
is not yet clear whether the coincidence of the lines with and without electric-current in the fluid 
occurs only for the considered model with a very narrow gap and perfectly-conducting cylinders 
or not. Note that in narrow gaps the radial profiles of the azimuthal fields between the cylinders 
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both are almost uniform. Indeed, test calculations also provided instability even for fields with 
uniform D cp in the same domain of Reynolds number and Hartmann number. One could believe 
that for Pm ^ 1 the super-rotation becomes unstable under the mere presence of any toroidal field 
but for Pm = 1 the dissipation processes prevent the excitation of this rather slow (see below) 
instability. 

B. Pattern migration 

In general both modes with m — ±1 are simultaneously excited. If the pinch rotates rigidly then 
both modes have exactly the same amplitudes and form a standing wave. The instability pattern 
looks azimuthally dipolar with the drift direction depending on the magnetic Prandtl number. Note 
also the nearly circular geometry of the resulting cells in the R-z -plane (Fig. [5]). For a certain 
Pm between 0.1 and 1 the azimuthal migration disappears and the entire pattern will rest in the 
laboratory. For nonuniform rotation laws with finite shear one of the modes m = 1 or m = - lis 
preferred and the instability pattern approaches a spiral. 



FIG. 5. The isolines of the bn in a narrow gap for Pm = 1 (left) and Pm = 0.1 (right) for rigid rotation. 
Both modes with m = ±1 are excited with the same amplitude forming a standing wave. The pattern with 
Pm = 1 migrates in the rotation direction while it migrates opposite for Pm = 0.1 . r- m = 0.95, // = I. 
Re = 111, Ha = 3440, perfect-conducting boundaries. 

Figure [6] shows the behavior of the azimuthal migration of the nonaxisymmetric vortices as 
more diverse. It is striking that for large Pm and/or for sub-rotation often |cudr| — l2 out . From Fig. 
|6] we also find that the nonaxisymmetric instability pattern for sub-rotation nearly corotates with 
the outer cylinder (as it is observed for AMRI 4,17 ) for all magnetic Prandtl numbers. For Pm > 1 
the drift frequencies for super-rotation are also negative so that their magnetic pattern azimuthally 
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migrates in die direction of the rotation. 



0 2500 5000 7500 10000 

HARTMANN NUMBER 


0 2500 5000 7500 10000 

HARTMANN NUMBER 



FIG. 6. The same as in Fig. [3]but for the normalized drift frequency u;,-] r / J? out (see Eq. (IT2l ) ) in a narrow gap 
(r in = 0.95) for Pm = 10 (left), Pm = 1 (middle) and Pm = 0.1 (right). r ui = 0.95, perfect-conducting 
boundaries. 


For smaller magnetic Prandtl numbers, however, for rigid rotation and for super-rotation the 
pattern counterrotates. The azimuthal migration of the linear solutions directly reflects the actual 
value of the magnetic Prandtl number. This is also true for a rigidly rotating pinch. After the results 
plotted in Fig. [6]its pattern corotates with the outer cylinder for small Pm and it counterrotates with 
the outer cylinder for large Pm. Hence, numerical simulations with magnetic Prandtl number unity 
for flows with vanishing or positive shear may easily lead to results which are not representative 
for the solutions with smaller Pm. On the other hand, rotation laws with negative shear do not 
show that sensitivity to the Pm-value (see also Fig. [7] below). 


V. WIDE GAP 
A. Instability map 

For a wide gap with r in = 0.5 and also in advance to a possible laboratory experiment Figs. [7] 
give the eigenvalues for marginal instability, the wave numbers and the drift rates for a fluid with 
the magnetic Prandtl number of Pm = 10 -5 (liquid sodium). The characteristic Hartmann number 
for conducting boundaries and for resting cylinders is HaTay = 35.3 (see Fig. [D. From now on 
Reynolds numbers and drift rates are related to the rotation rate of the outer cylinder. 

Figure [7] (left panel) shows the curves of marginal instability for many rotational laws with // 
between 1 and 128. The subcritical excitation of the current-driven instability by super-rotation is 
much stronger than for rigid rotation (dashed curve). If the Reynolds number is formed with the 
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FIG. 7. Left panel: instability map for super-rotation in a wide gap, the lines are marked with \x. The 
dashed curve belongs to rigid rotation while for reasons of comparison the dotted curve shows the result for 
// = 0.25 (Rayleigh limit). Middle panel: drift rates cu ( ] r /i? ou t. Right panel: the vertical wave numbers are 
nearly uniform. r in = 0.5, Pm = 10 -5 , perfect-conducting cylinders. 


outer rotation rate, i.e. Re out = /iRe, then the curves for strong shear are converging. Hence, a 
minimum Hartmann number of order 30 exists which cannot further be reduced by steeper rotation 
laws. It is also clear that for Re out > 200 the differential rotation starts simply to suppress the 
nonaxisymmetric instability. The middle panel of the Fig. [7] shows that the azimuthal drift of 
the instability pattern has the same (positive) sign for rigid rotation and super-rotation and the 
opposite (negative) sign for sub-rotation. The idea is thus supported that the phenomenon of 
subcritical excitation, i.e. dRe/dHa < 0, for slow rotation with vanishing or positive shear is 
a common double-diffusive phenomenon which disappears for Pm = 1. On the other hand the 
known lowering of the critical Hartmann number and the negative drift values by sub-rotation 
exists for all Pm (see Fig. [2]). Note that for ccdr/-O olll = —1 the magnetic pattern strictly corotates 
with the outer cylinder. In contrast, the magnetic pattern for all super-rotation laws including 
uniform rotation migrates opposite to the sense of rotation corresponding to the behavior of the 
drift rates in the narrow gap. For small Pm by choice of the rotation rate of the outer cylinder one 
can obtain all sorts of migration of the magnetic pattern between corotation and counterrotation. 

The axial wave numbers k are plotted in the right panel of Fig. [T] They all show similar values. 
Let Sz be the characteristic vertical scale of a cell. Then from the definitions follows 


5z 

15 


-Rin 

D 


(13) 


so that 8z/D ~ n/k forms the relation between the cell size and the wave number for r in = 0.5. 
Hence, the data with k ~ tt lead to 8z ~ D what means that the Tayler cells approximately form 
a circle in the R-z plane. In this general formulation the result does not depend on the gap width. 
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B. Growth rates 


The subcritical excitation of the TI for rotation laws with positive shear only exists for suffi¬ 
ciently slow rotation (Mm 1) so that the instability only grows very slowly as Fig. [8| demon¬ 
strates for super-rotation. The growth rates are normalized with the outer rotation rate. Then 
the maximum growth rates always occur for the same Reynolds number. For the upper curve 
of the plot one finds that ce gr ~ 0.03i2 O ut for Re ~ 130 so that the exponential growth time is 
Tgr ~ 10f? 2 ut in seconds when R oxlt is measured in cm (it is v = 7 ■ 10~ 3 cm 2 /s for liquid sodium). 
If the expression 


d2 

w gr = r- out 


fJ-oPV 


( 14 ) 


of the growth rate is adopted (with T as a dimensionless numerical factor which only depends 
on the gap width and the magnetic Prandtl number), which has been derived for the nonrotating 
pinch and which has been experimentally realized~54& 5 then for the strongest super-rotation the 
Fig. [8] gives T ~ 10~ 3 which very well fits the theoretical results for the resting container. This 
value certainly increases for wider gaps but the exponential growth time of the instability for slow 
super-rotation will hardly be shorter than that for the resting pinch. 



FIG. 8. Growth rates normalized with i? out for Ha = 33 with // = 4, 8, and 128. r m = 0.5, Pm = 10 5 , 
perfect-conducting cylinders. 
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VI. SUMMARY AND OUTLOOK 


A rotating pinch with a homogeneous axial electric-current has been considered where the 
cylindric bounding walls rotate with /x > 1 , i.e. the outer cylinder rotates with the same rotation 
rate or rotates faster than the inner cylinder. A linear perturbation theory fixes the critical Hart¬ 
mann numbers for which the system becomes marginally unstable. The surprising result is the 
occurrence of a double-diffusive instability which only exists for v / 77 . For slow rotation with 
l7 ou t f?A the excitation of the nonaxisymmetric perturbations becomes subcritical, i.e. the crit¬ 
ical Hartmann number for rotation is smaller than without rotation. The effect is rather weak for 
rigid rotation but it is remarkably strong for super-rotation. It is numerically shown that for steeper 
and steeper rotation laws the series of minimum Hartmann numbers converge to a total minimum 
HaMin for n —)■ 00 (approaching resting inner cylinders). The resulting normalized lowering of Ha 

Ha^ay HaMin 

X= -„- ( 15 ) 

HaTay 

depends on the magnetic Prandtl number Pm. It vanishes for Pm = 1 and takes similar values for 
very large and for very small Pm. Relaxations of the the critical Hartmann number of order 20% 
(rin = 0.5) and 80% (r in = 0.95) exist for Pm 7 ^ 1 (Fig. left panel). Note also that the excess 
m grows for decreasing gap width. 

It is also shown for a narrow gap and for small magnetic Prandtl number that in the area of 
the instability map where the subcritical excitation exists for super-rotation the form of the lines 
of marginal instability only weakly depends on the radial profile of the azimuthal magnetic field. 
For finite Reynolds number the lines of marginal instability in Fig. [4] for the two different profiles 
© and © are very close together. In this domain of the map the radial distribution of the axial 
electric-current seems to be unimportant. Obviously, the energy provided by the differential rota¬ 
tion is large enough to maintain the instability while the magnetic field is only needed as a catalyst. 
This is a numerical finding which implies that among all other radial profiles also the profile with 
B r oc 1 /R (no current within the fluid) leads to a nonaxisymmetric instability for super-rotation. 
Again, however, this phenomenon disappears for Pm = 1 revealing its double-diffusive character. 

The instability for Pm 1 and for Pm -C 1 differs in another respect. For sub-rotation we 
always find that the pattern migrates for all Pm in positive direction of the azimuthal coordinate </>. 
For solid-body rotation and for super-rotation there is, however, a strong influence of the magnetic 
Prandtl number on the azimuthal migration of the perturbation patterns. For fi > 1 the pattern 
counterrotates for small Pm while it corotates for Pm > 1. The Pm-dependence of the drift 
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FIG. 9. The maximal reduction of the Hartmann number (left) and the related drift frequencies cu dr /l7 out 
(right) as function of the magnetic Prandtl number. For Pm 1 and for Pm <C 1 the subcritical excitation 
measured by (fl5l) is very similar but the signs of the azimuthal migration of the instability pattern differ. 
The instability patterns drift in the rotational direction (corotation) or opposite (counterrotation). Note the 
strong influence of the gap width on the Hartmann number reduction. The curves are marked with r; n , 
H =—> oo, perfect-conducting boundaries. 

frequency cu dr disappears for small and for large magnetic Prandtl numbers (Fig. [9] right panel). 
The latter finding is also true for a pinch with solid-body rotation. 

A final question is whether the described phenomena can be realized in the laboratory. The 
answer is yes and there is an interesting variety of possibilities.The following estimations are 
derived for liquid sodium with Pm = 10 -5 as the fluid conductor. In Fig. [7] for a container with 
r in = 0.5 the characteristic Hartmann number for resting containers is 35.3 while the smallest 
Hartmann number for (strong) super-rotation is 30.4. The necessary current in the gap to produce 
a certain Hartmann number Ha is 

r c /(l-r in )(l + r in ) 2 ,- u 

-'fluid = 5w - 3 - y/tM)PVV Ha (16) 

y An 

(Rudiger et al. 2013—) which does not depend on the physical size of the container. For liquid 
sodium it is y/noP ur l — 8-2 G-cm (for gallium ~ 26 G-cm). To reach Ha ~ 30 one needs an 
electric-current of 3.6 kAmp flowing through the sodium (11 kAmp for gallium). 

After Fig. 0] experiments are more interesting for narrow gaps with r in = 0.95. A character¬ 
istic minimum Hartmann number is 0(1000) for super-rotation which needs about 20 kAmp for 
its generation with liquid sodium. The Tayler instability with resting cylinders would, however, 
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require about 45 kAmp for its realization which seems to be much too high. Only the interplay 
with differential rotation allows this instability to observe. A particular challenge for such experi¬ 
ments is that already for a slight increase of this Hartmann number the electric-current through the 
sodium can be replaced by an axial current inside the inner cylinder. Both resulting instabilities 
are very similar with respect to wave number and azimuthal drift. Such an experiment would eas¬ 
ily demonstrate the surprisingly close relation of AMRI and TI under the presence of differential 
rotation. 

As the instability phenomena presented in this paper all belong to the class of small magnetic 
Mach numbers the technical realization of the Reynolds numbers required by the lines in Figs. |4] 
and|7]is also no problem. Note, however, that the needed rotation rates of the cylinders do indeed 
depend on the physical size of the experiment. 
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Appendix A: A local approximation 

The described effects have mainly been calculated with boundary conditions for perfect¬ 
conducting cylinders. After Fig. [2] the negative slope dRe/dHa < 0 of the lines of marginal 
instability for slow rotation also exists for models with vacuum conditions but with reduced effi¬ 
ciency. It is thus worthwhile to discuss in plane geometry (for narrow gaps) the result 

pc * 2 _ 1 ((1 + Ha* 2 m* 2 ) 2 - 4Ha* 4 m* 2 )(l + Ha* 2 m* 2 ) 2 

“ 4 Ha* 4 Ro 2 m* 2 - (Ro + 1)((1 + Ha* 2 m* 2 ) 2 - 4Ha* 4 m* 2 ) ’ 

of a local approximation for Pm — > 0 where Re*, Ha* and m* represent the slightly modified 
Reynolds number, Hartmann number and azimuthal wave number—. The Rossby number Ro = 
(l/2)dlogi?/dlogi? represents the differential rotation, it is positive for super-rotation and negative 
for sub-rotation. In contrast to all other quantities the Rossby number enters the Eq. (IA1I) with odd 
and even powers. In (IA1I) the Rossby number Ro must be independent of the radius R which only 
allows to consider Taylor-Couette flows of 0.2 < fi < 2 (for r in = 0.5). The latter (super-Rotation 
law corresponds to Ro ~ 1 in a very good approximation. 
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We shall show that all curves with m* > 1 in a Ha*-Re* plane for slow rotation show a subcrit- 
ical behavior compared with the critical eigenvalue 11 a j, iy = 1/ y/m*( 2 — m*) for Re* = 0, i.e. 
Ha*(Re*) < Ha^ y . To this end the function Z — (1 + Ha* 2 m* 2 ) 2 — 4Ha* 4 m* 2 is defined so that 
(1A1I) yields 

4Re* 2 Ha* 4 Ro 2 m* 2 

~~ (1 + Ha* 2 m* 2 ) 2 + 4(Ro + l)Re* 2 ' 

For rigid rotation (Ro = 0) only the solution Z = 0 exists which does not reflect the rotational 
influence as given by Fig. [1] Obviously, the function Z(Re*) only vanishes for Re* = 0 and 
it is positive-definite for finite Re* if - as we shall assume - Ro > —1. The above mentioned 
eigenvalue Ha^ ay forms the solution of Z — 0. It only exists for m* < 2. The solution of Z = 5 
with S > 0 can thus be written as Ha* 2 = Ha^ 2 y + £ with unknown sign of e which, without loss 
of generality, can be assumed as small against Ha.y ay . Hence, from the definition of the function 
Z follows 


5 

4 m* ’ 


(A3) 


so that always £ < 0. For Ro > —1 it is thus Ha*(Re) < Ha^ for negative and positive shear, 
i.e. the excitation of the TI becomes always subcritical by the action of any differential rotation. It 
is possible to demonstrate that this result does not change without the restriction to small e. 

The fact that Z = 0 requires Re* = 0 has the consequence that Ha* = Ha^. lv does not appear as 
a solution of (IA1I) for finite Reynolds numbers. Hence, the curves of marginal instability always 
remain in the subcritical domain with Ha* (Re *) < Ha^ ay and never reach Hartmann numbers 
larger than Ha^ ay . The typical suppression of the magnetic instabilities by fast rotation (see Fig. 
[7]) is thus not reflected by the local relation (IA1I) for inductionless fluids. 
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